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Abstract 



We discuss the physics underlying event-by-event fluctuations in rel- 
ativistic heavy ion collisions. We will argue that the fluctuations of 
the ratio of positively over negatively charged particles may serve as a 
unique signature for the Quark Gluon Plasma. 

> 

: 1 INTRODUCTION 

O 

Any physical quantity measured in experiment is subject to fluctuations. In 
psj ! general, these fluctuations depend on the properties of the system under study 

(in the case at hand, on the properties of a fireball created in a heavy ion 
collision) and may contain important information about the system. 

The original motivation for event-by-event (E-by-E ) studies in ultra rei- 
ativistic heavy ion collisions has been to find indications for distinct event 
classes. In particular it was hoped that one would find events which would 
carry the signature of the Quark Gluon Plasma. First pioneering experiments 
in this direction have been carried out by the NA49 collaboration Q. They 
analyzed the E-by-E fluctuations of the mean transverse momentum as well as 
the kaon to pion ratio. Both observables, however, did not show any indication 
for two or more distinct event classes. Moreover, the observed fluctuations in 
both cases were consistent with pure statistical fluctuations. 

On the theoretical side, the subject of E-by-E fluctuations has recently 
gained considerable interest. Several methods to distinguish between statistical 
and dynamic fluctuations have been devised 0, [|. Furthermore the influence 
of hadronic resonances and possible phase transitions has been investigated 
[f§, |f| H, [7|, §. All these theoretical considerations assume that the observed 
fluctuations will be Gaussian and thus the physics information will be in the 
width of the Gaussian, which is controlled by 2-particle correlations ||. 



2 EVENT-BY-EVENT FLUCTUATIONS 



Fluctuations have contributions of different nature. Besides the statistical 
fluctuations due to a finite number of particles in case of heavy ion collisions, 
there are also fluctuations of the volume. Both these fluctuations are trivial 
and add to the dynamical fluctuations which carry the real information about 
the properties of the system. The dynamical fluctuations are controlled by 
the appropriate susceptibilities, which are the second derivative of the free 
energy with respect to the appropriate conjugate variable. For example the 
fluctuations of the charge are given by 

m 2 ) = -T0 = -VT XQ (1) 

where [iq is the charge chemical potential, T the temperature and V the volume 
of the system. \q is the charge susceptibility. It is interesting to note that 
the very same susceptibility also controls the response of the system to an 
external electric field. The properties of a macroscopic system are studied by 
investigating its response, i.e. its susceptibility, to external forces. This is of 
course impossible for the mesoscopic systems created in a heavy ion collision. 
However, the same property can also be accessed via fluctuations. 

Since we will mostly concentrate on the charge susceptibility, let us point 
out that they are directly related to the electric mass, which is given by the 



zero momentum limit of the static current-current correlation function |T_0 

n 00 (^ = o,£^o). 

(5Q) 2 ) = -VTU 00 (uj = 0,k^0) (2) 



This object has been, evaluated in Lattice QCD [11, [L2|] as well as in effec- 



tive hadronic models (see below and ]TD|, |Tj| |T3[). Note also that this is the 
same current-current correlation function, which controls dilepton and photon 



production in heavy ion collisions JTfl, 17 . 



2.1 Fluctuations of particle ratios 

In order to avoid volume fluctuations one needs to study observables which are 
independent of the volume of the system. Among others the ratio of particle 
multiplicities will have this property. This is certainly true if one looks at 
similar particles such as 7r + and ir~ , where the freeze out volumes are expected 
to be the same. Some residual volume fluctuations may be present if one 
considers ratios of particles with different quantum numbers such as the K/tt 



ratio, but they still should be small. Let us define the particle ratio R\ 2 of two 
particle species N\ and N 2 

The fluctuations of this ratio are then given by [Q, §, [| 

(4) 



(5R 12 ) 2 ( ((SN,) 2 ) ((5N 2 ) 2 ) (5N 1 5N 2 y 
(R12) 2 ^ V (Nif (N 2 ) 2 WW. 



The last term in Eq. (Q) takes into account correlations between the particles of 
type 1 and type 2. This term will be important if both particle types originate 
from the decay of one and the same resonance. For example, in case of the 
7T + /7r _ ratio, the p$, w etc. contribute to these correlations. Also this term is 
responsible to cancel out all volume fluctuations ||. 

Let us note that the effect of the correlations introduced by the resonances 
should be most visible when (Ni) ~ (-^2)- On the other hand, when (N 2 ) ^> 
(iVi), as in the K/tt ratio, the fluctuations are dominated by the less abundant 
particle type and the resonances feeding into it. The correlations are then very 
hard to extract. In || it was shown that in case of the K/ir-rdXio resonances 
and quantum statistics give rise to deviations from the statistical value of at 
most 2 %, in agreement with experiment JlJ. 

As pointed out in [|J the measurement of particle ratio fluctuations can 
provide important information about the abundance of resonances at chemical 
freeze out, and thus provides a crucial test for the picture emerging from 
the systematics of single particle yields ||18|| . In particular the fluctuations 
of the 7r + /7r~-ratio should be reduced by about 30 % as compared to pure 
statistics due to the presence of hadronic resonances with decay channels into 
a 7r + -7r~-pair at chemical freeze out. About 50 % of the correlations originate 
from the decay of the po an d the u mesons. Thus the fluctuations provide a 
complementary measurement to the dileptons. 



3 Fluctuations and correlations 

As already pointed out in the previous section, the fluctuations are sensitive 
to correlations among the particles of concern. In case of charge fluctuations 
which we will discuss below, the relation between the charge fluctuations and 
the 2-particle correlation function C 2 for a system of unit charged particles is 
given by H 

((5Q) 2 ) = C 2 (+, +) + C 2 (- -) - 2C 2 (+, -) + Pl (+) + pi(-) (5) 



where the correlation function C 2 is defined as 

d 2 N 



Co(±,±) 



dyi,±dy 2i ±- -, 

Ay dyi t ±dy 2 ,± 



dN dN 
dyi,± dy 2 .± 



The single particle density pi(±) is defined as 



Pi(±) 



, dN 
dy±-7— 

Ay dy± 



(6) 



(7) 



Here Ay denotes the range of acceptance e.g. in rapidity, over which the 
particles are measured. This expression can easily be extended to a system of 
particles with charge states different from unity . 

Thus, an alternative way to access charge fluctuations is the measurement 
of the one and two particle densities. However in this case the effect of volume 
fluctuations still needs to be removed. This can be achieved by defining appro- 
priate modifications of the correlation functions as discussed in [2T|. Also the 
so-called balance functions introduced in ref. |22 can be viewed as another 
way of analyzing these correlation functions and a relation between the ratio 
fluctuations discussed here has been made in [23|. 



4 Charge fluctuations 

Measuring the charge fluctuations or more precisely the charge fluctuations per 
unit degree of freedom of the system created in a heavy ion collision would tell 
us immediately if we have created a system of quarks and gluons |24j (see also 
p5fl). The point is that in the QGP phase, the unit of charge is 1/3 while in 



the hadronic phase, the unit of charge is 1. The net charge, of course does not 
depend on such subtleties, but the fluctuation in the net charge depends on 
the squares of the charges and hence are strongly dependent on which phase it 
originates from. However, as discussed in the previous section, measuring the 
charge fluctuation itself is plagued by systematic uncertainties such as volume 
fluctuations, which can be avoided if one considers ratio fluctuations. The task 
is then to find a suitable ratio whose fluctuation is easy to measure and simply 
related to the net charge fluctuation. 

As shown in the fluctuation of the ratio R = N + /N_ serves this purpose 



and the observable to investigate is 



D ee (AT ch ) (W? 2 ) = 4 (AT*) (5F 2 ) = 4^1 (8) 



which provides a measure of the charge fluctuations per unit entropy. 



In order to see how this observable differs between a hadronic system and 
a QGP let us compare the value for D in a pion gas and in a simple model of 
free quarks and gluons. Ignoring small corrections due to quantum statistics 
a simple calculation gives for a pion gas 

Dji—gas ~ 4 . (9) 

In the presence of resonances, this value gets reduced by about 30 % due to 
the correlations introduced by the resonances, as discussed in section (|2.1| ). 
For a thermal system of free quarks and gluons we have in the absence of 



correlations and ignoring small correction due to quantum statistics |24| 



A3GP-3/4. (10) 

Actually the charge fluctuations ((5Q) 2 ) have been evaluated in lattice 
QCD along with the entropy density |TT[. Using these results one finds 



y Lattice-QCD 



1.5 fir 



where the uncertainty results from the uncertainty of relating the entropy to 
the number of charged particles in the final state. Actually the most recent lat- 



tice result |L2| for the charge fluctuations, which was obtained in the quenched 
approximation, is somewhat lower then the result of [O]. 

But even using the larger value of D = 1.5 for the Quark Gluon Plasma, 
there is still a factor of 2 difference between a hadronic gas and the QGP, 
which should be measurable in experiment. 

The key question of course is, how can these reduced fluctuations survive 



the hadronic phase. This has been addressed in |24| and in more detail in p6 
The essential point why this signal should survive is that charge in conserved 
which leads to a very slow relaxation of the initial fluctuations. 

Finally one needs to take into account that the total charge of the system 
in conserved and thus does not fluctuate. This has been addressed in and 
a corrected observable D has been derived. 

In Fig. ([I]) the importance of these corrections, in particular the charge 
conservation correction is demonstrated. There we compare the uncorrected 
observable D with the corrected observable D as a function of the width of 
the rapidity window based on a URQMD [^8] simulation. The effect of the 



charge conservation is clearly visible. With increasing rapidity window the 
charge fluctuation are suppressed. Once the charge-conservation corrections 
are applied, the value for D remains constant at the value for a hadron gas of 
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Figure 1: Charge fluctuations as a function of the size of the rapidity window. 



D ~ 3. This is to be expected for the URQMD model, which is of hadronic 
nature and does not have quark- gluon degrees of freedom. 

Also for small Ay the value is D ~ 4 before it drops to D ~ 3 for Ay > 1.5. 
This effect, which was already predicted in 0, is simply due to the fact that 
the correlation introduced by the resonances gets lost if the acceptance window 
becomes too small. 



Recently an interesting comparison |29[ between predictions for D by dif- 
ferent event generators has been made. The authors find that indeed a parton 
cascade model leads to small values for D provided hadronic rescattering is 



turned off. This is in line with the predictions of P^| . However, the authors 
also report substantial corrections due to hadronic rescattering, which appear 



quite large given the arguments presented in ]24|, [26j . But none of the event 
generator gave result large than D ~ 3. 

5 Effects of particle interactions 

So far we have discussed the fluctuations of free, noninteracting particles only. 
Of course, by considering resonances as well, some of the interactions are 
taken into account [14, 15]. To address the effect of interactions we have 



calculated the electromagnetic current-current correlator II 00 as well as the 
entropy to the first order of the interaction for a system of pions interacting via 
p-meson exchange. For simplicity, we have assumed that the p-mass is infinite, 
which leads to point-like interaction among the pions, and is goverened by the 



following effective Lagrangian 



4 ^ 



2i ^+i((w a -^)+i^r 



2mj; 







4o (^)* - < 

The 2-loop diagrams for the corrections to n o are shown in fig I 
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Figure 2: Two loop diagrams contributing to n o to order g 2 . 

The ratio H 00 (T)/S(T) for the entropy S(T) and the self energy n 00 (T) is 
shown in fig [5| 

The effect of the interactions is to increase the charge fluctuations per 
entropy. At temperature of ~ 140 MeV the corrections are quite large, in line 
with the findings of U|, where chiral perturbations theory has been used. 
However, in the same paper JTH] the effect the charge fluctuations have also 
been evaluated in a virial expansion, leading to much smaller (but still positive) 
corrections. Therefore, it appears that at least dynamical p mesons need to be 
taken into account before any quantitative conclusion can be made. 



6 Conclusions 

We have discussed event-by-event fluctuations in heavy ion collisions. These 
fluctuations may provide useful information about the properties of the matter 
created in these collisions, as long as the 'trivial' volume fluctuations, inherent 
to heavy ion collisions, can be removed. We have argued that the fluctuations 
of particle ratios is not affected by volume fluctuations. 
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Figure 3: Ratio of charge fluctuations over entropy for free particles (full line) 
and with interactions (dashed line). 



In particular the fluctuations of the ratio of positively over negatively 
charged particles measures the charge fluctuation per degree of freedom. Due 
to the fractional charge of the quarks, these are smaller in a QGP than in a 
hadronic system. 

A measurement of our observable D ~ 1 would provide strong evidence 
for the existence of a QGP in these collisions. A measurement of D ~ 3 on 
the other hand does not rule out the creation of a QGP. There are a number 



of caveats (see f24|), which could destroy the signal, such as unexpected large 



rapidity shifts during hadronization. 

At this workshop first preliminary results for the results observable D have 
been reported by the STAR, CERES and NA49 collaboration |U|. They all 
find a value which is consistent with D ~ 4. This finding is somewhat surpris- 
ing since not even the effect of resonances seem to be visible in the fluctuations. 
At the same time STAR reports that the measured number of K* mesons is 



well in line with the expectations of the thermal model [31, 19|. So were are 



the correlations? Could it be that the interactions among the resonance are 
indeed so strong that they entirely compensate the reduction of D due to res- 
onances. Also, as already mentioned, in the study presented in |29| none of 
the event generators finds a value of D > 3 for rapidity windows Ay > 2. This 
certainly needs further investigation. 

Finally let us note that fluctuations of the baryon number in principle can 
also be utilized, since in the QGP quarks carry fractional baryon number p5| . 



This, however, would require the measurement of neutrons on an event by 



event basis. 
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